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ABSTRACT 

'VWVW\A* 

A theoretical  model  is  developed  to  describe  the  synergic  effects  of 

» 

analyte  vaporization  and  diffusion  rates  on  the  spatial  distribution  of 
analyte  atoms  in  the  vapor  cloud  surroundirig  an  individual  solute  particle 
in  a laminar  flame  or  plasma.  A mathematical  convolution  of  the  vapori- 
zation and  diffusion  processes  is  used  to  describe  the  analyte  spatial 
distribution.  Unfortunately,  the  resulting  convolution  integral  cannot 
be  solved  in  closed  form  for  the  experimentally  observed  vaporization 
function,  and  several  simplifying  approximations  to  the  vaporization  func- 
tion are  compared  and  shown  to  be  useful.  With  this  new  model,  the 
analytical  consequences  of  changes  in  both  diffusion  coefficient  and 
vaporization  rate  are  examined,  and  the  resulting  changes  in  spatial  dis- 
tribution of  atoms  about  a typical  analyte  particle  are  compared. 


. 

An  important  endeavor  in  the  improvement  of  all  atomic  .spectromei ric 
methods  of  analysis  is  the  study  of  spatial  atom  distributions  in  tin- 
analytical  source.  A quantitative  knowledge  of  the  physical  and  chemical 
factors  that  determine  a particular  free  atom  distribution  is  Important 
in  understanding  the  influence  of  many  interferents  and  is  an  essential 
part  of  the  effort  to  achieve  absolute  flame  or  plasma  atomic  spectrometric 
measurements  (1). 

The  distribution  ot  atoms  in  an  analytical  flame  or  plasma  strongly 
affects  the  optimal  experimental  configuration  to  be  employed  in  atomic 
spectrometry.  For  example,  in  flame  spectrometry,  an  optimal  flame  height 
exists  where  most  of  the  analyte  atoms  have  been  released  but  where  they 
have  not  yet  fully  equilibrated  through  ionization  or  recombination 
reactions.  In  addition,  the  horizontal  distribution  of  atoms  across  a 
flame  or  plasma  determines  the  proper  optical  aperture  for  measurement  of 
the  atoms,  be  it  by  absorption,  emission,  or  fluorescence.  For  instance, 
a narrow  atomic  spatial  distribution  in  an  atomic  absorption  experiment 
will  cause  greater  non-linearity  of  working  curves  due  to  source  radiation 
entering  the  detector  without  having  passed  through  the  sample  atoms.  In 
an  atomic  fluorescence  experiment  as  well,  a narrow  atomic  vapor  cloud 
would  experience  a smaller  cross-section  of  the  exciting  radiation  than  a 
broader  cloud.  Lastly,  in  atomic  emission  spectrometry,  a high  concentra- 
tion of  analyte  atoms  within  a narrow  region  of  the  flame  or  plasma  center 
would  produce  different  self-absorption  characteristics  than  those  produced 
by  the  same  number  of  atoms  distributed  over  the  entire  source. 

Conveniently,  the  same  processes  affect  the  distribution  of  atoms  in 
both  flames  and  plasmas,  so  that  studies  designed  to  explore  these  phenomena 
in  one  source  will  yield  results  applicable  to  the  other.  In  both  kinds  of 
source,  final  spatial  variations  in  atom  concentration  are  governed  by  the 
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initial  distribution  of  aerosol  droplets,  by  the  location  of  the  dried 
(and  perhaps  fragmented)  solute  particles  which  result  from  dosolvation  of 
the  original  aerosol,  and  by  the  spreading  of  atomic  vapor  produced  upon 
vaporization  of  the  solute  particles.  In  turn,  localized  variations  in 
the  number  density  of  droplets,  particles,  and  atoms  are  controlled  by 

two  distinguishable  forces:  convective  and  diffusive.  Convective  forces  j 

arise  from  motion  of  the  flame  or  plasma  gases  about  the  droplets  or  par- 
ticles while  diffusive  forces  arise  from  concentration  gradients  and  serve 
to  spatially  spread  the  droplets,  particles,  and  atoms  apart  from  each 
other.  Clearly,  prediction  of  the  final  atom  concentration  in  an  analyti- 
cal source  will  require  a complete  understanding  of  the  effect  of  both 
these  forces  on  the  sample  components  during  the  atomization  process. 

In  flame  atomic  spectrometry,  it  has  been  established  that  contami- 
nants or  matrix  components  in  nebulized  solutions  can  cause  changes  in 
the  lateral  distribution  of  analyte  atoms  in  the  measurement  region  (2-4). 

Clearly,  this  kind  of  interference  can  be  quite  troublesome  in  practical 
flame  spectrometry.  For  example,  if  a sample  and  standard  solution  create 
different  lateral  atom  distributions  and  the  spectrometer's  optical  system 
is  optimized  for  the  distribution  of  the  standard,  the  sample  measurement 
will  produce  either  lower  or  higher  sensitivity,  depending  on  whether  the 
sample  atoms  are  more  or  less  dispersed  than  those  of  the  standard. 

L’-'ov  and  co-workers  (5-7)  have  attributed  these  lateral  distribution 
anomalies  to  differences  in  the  size  and  density  of  the  aerosols  generated 
from  the  two  solutions,  and  a consequent  change  in  the  vertical  spread  of 
the  aerosols  as  they  emerge  from  the  burner  top.  Tn  contrast.  West,  Fassel, 
and  Kniseley  (8)  ascribed  the  lateral  df.str.lbul Ion  vartal Ions  to  an  inter- 
fered -produced  change  In  the  maximum  vapor  I y.al  * >>n  region  of  1 he  analyte. 

In  practice,  either  mechanism  for  the  interference  seems  viable  and  In 
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west  situations  the-  two  probably  act  in  concert  to  determine  the  observed 
atom  distributions.  in  order  to  understand  how  the  lateral  spread  of 
aerosol  particles  and  the  vaporization  rate  of  those  particles  interact  to 
produce,  an  observed  lateral  distribution  of  atomic  vapor,  and  indeed  to 
understand  the  detailed  microscopic  processes  leading  to  atom  formation  and 
distribution,  it  is  most  meaningful  to  examine  the  spatial  distribution  of 
atoms  created  by  an  individual  aerosol  particle  in  a laminar  analytical 
flame . 

Accordingly,  the  objective  of  this  paper  is  to  develop  an  accurate 
model  of  tins  spatial  distribution  of  atoms  generated  by  volatilization  of 
an  individual  solute  particle  in  a laminar  analytical  flame.  This  model 
portrays  the  combined  effects  of  vaporization  and  diffusion  as  a convolu- 
tion if  the  two  processes  and  neglects  the  influence  of  vapor-phase  reac- 
tions on  the  resulting  spatial  atomic  distributions.  A rigorous  numerical, 
solution  to  this  convolution  will  be  presented  and  evaluated.  In  addition, 
approximate  vaporization  rate  functions  will  be  compared  which  allow  analyti- 
cal solutions  to  the  diffusion-vaporization  convolution,  and  which  yield 
more  readily  utilizable  mathematical  expressions.  In  a paper  to  follow, 
experimental  measurements  will  be  used  to  verify  the  validity  of  the  work 
and  model  presented  herein. 

THEORETICAL  MODEL 

WWlA/'/MAAn/Vm 

As  a solution  droplet  travels  through  an  analytical  flame,  the  heat 
transferred  to  it  from  the  flame  causes  the  solvent  to  evaporate  (9).  Once 
desol  vat  ion  is  complete,  the  surface  temperature  of  the  dried  aerosol 
particle  will  rise  and  soon  reach  a level  sufficient  to  cause  noticeable 
volatilization  of  analyte  atoms.  Because  these  events  arc  kinetically 
coni  rolled,  droplets  of  n part  icular  size  and  which  follow  a particular 


path  through  the  flame  will  require  a reproducible  residence  time  in  the 
flame  before  analyte  vaporization  begins  (10).  Correspondingly,  because 
droplets  and  their  vapor  products  travel  vertically  with  the  flame  gases, 
any  variation  in  the  time  required  for  dcsolvation  will  merely  produce  n 
change  in  the  flame  height  where  vaporization  begins.  Other  than  this 
predictable  vertical  shift,  changes  in  droplet  residence  time  will  have 
little  effect  on  the  spatial  distribution  of  atoms  created  by  an  individual 
aerosol  droplet  and  will  be  ignored  here.  Accordingly,  in  the  present 
model,  the  time  when  volatilization  of  analyte  atoms  first  begins  will  be 
assigned  the  time  coordinate  reference  _t  = 0;  any  displacement  of  this 
temporal  reference  can  then  be  easily  incorporated  into  further  treatments. 

As  soon  as  an  solution  droplet  is  introduced  into  a laminar  flame,  the 
rapidly  rising  flame  gases  act  through  viscous  forces  to  accelerate  it; 
by  the  time  the  droplet  is  completely  desolvated,  its  velocity  is  experi- 
mentally indistinguishable  from  that  of  the  flame  (10-12).  In  turn,  because 
a dried,  vaporizing  solute  particle  travels  essentially  at  the  flame  gas 
velocity,  it  can  be  assumed  to  be  surrounded  symmetrically  by  its  vapori- 
zation products.  Consequently,  a spherical  spatial  coordinate  system  can 
conveniently  be  adopted,  whose  origin  (r_  = 0)  lies  at  the  center  of  mass 

t. 

of  the  vapor  cloud.  Of  course,  because  the  cloud  travels  at  or  near  the 
flame  velocity,  this  origin  is  not  fixed.  However,  by  defining  such  a 
moving  reference,  the  spherically  symmetrical  cloud  can  be  described  in 
terms  of  a single  spatial  variable  jr. 

With  the  aid  of  this  spatial  and  temporal  coordinate  system,  the 
distribution  of  atoms  about  a volatilizing  particle  can  be  described. 

Analyte  atoms  released  during  the  vaporization  process  can  be  carried  away 
from  the  particle  by  cither  diffusion  or  convection.  Of  course,  in  a 
laminar  analytical  flame  such  an  most  often  employed  in  practical  flame 
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spectrometry,  the  predominant  mass  transfer  phenomenon  is  diffusion.  There- 
fore, the  transport  of  analyte  atoms  within  the  spherical  vapor  cloud  is 
governed  principally  by  Fick’s  laws  (13),  and  the  radially  and  temporally 
dependent  concentration  of  analyte  atoms  C(r, t) . will  be  given  by: 

* I a) 

In  equation  1,  _t  is  time  measured  from  the  onset  of  analyte  vaporization, 

D.  Is  the  diffusion  coefficient,  and  _r  is  the  radial  distance  from  the  center 
of  the  spherical  vapor  cloud. 

If  vaporization  were  instantaneous,  equation  1 could  be  solved  in  a 
rather  straightforward  manner.  To  do  so,  let  be  the  number  of  atoms 
vaporized;  then  at  any  time  after  vaporization, 


C(r , t) 


SOiDt)377 


-rJ/4Dt 

e 


(2) 


Unfortunately,  vaporization  occurs  on  the  same  time  scale  as  does  the  diffu- 
sive transport  of  atoms.  Thus,  each  successive  increment  of  analyte  vapor 
must  diffuse  into  a region  which  has  a concentration  of  analyte  atoms  already 
established  by  the  previously  vaporized  analyte.  Because  diffusive  forces 
are  controlled  by  concentration  gradients,  the  solution  of  Fick's  Laws  for 
this  problem  must  account  for  the  fact  that  analyte  vaporization  occurs 
during  the  diffusive  transport  process.  This  accounting  is  possible  through 
the  incorporation  of  a vaporization  function  4>(t)  into  equations  1 and  2. 

To  incorporate  a vaporization  function  into  equation  1,  a solution 
must  be  found  which  sums  the  effects  of  each  small  increment  of  analyte 
vapor  introduced.  The  solution,  accounting  for  the  synergic  effects  of 
diffusion  and  volatilization  in  creating  a spatial  distribution  of  atoms, 
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is  a convolution  and,  with  a generalized  analyte  vaporization  rate,  <l'(t), 
takes  the  form: 

, v 1 rtf  -rs/liD(t-t') 

c(r,t)  a J0  $(t,)  — fU’  (<) 

In  this  equation,  t/_  is  a delay  parameter  accounting  for  the  vaporizat  ion 

r 

from  the  beginning  of  vaporization  (where  t = 0)  and  summing  the  effects 
of  each  increment  up  to  the  last  time  to  be  considered  (where  t = t^) . To 
apply  equation  3 to  the  calculation  of  atomic  concentrations  at  different 
times  during  vaporization  requires  that  the  range  of  the  integral  (the  time 
range  under  consideration)  include  all  significant  analyte  vaporization; 
thus,  t^  must  be  set  equal  to  the  observation  time  _t . However,  when  post- 
vaporization times  are  of  interest,  only  the  vaporization  period  need  be 
included  in  the  range  of  the  integration,  so  that  tf  must  then  be  set  equal 
to  the  time  corresponding  to  the  completion  of  vaporization. 

Of  course,  for  the  above  treatment  to  accurately  describe  the  spatial 
distribution  of  analyte  atoms,  an  appropriate  vaporization  function  ^(t) 
must  be  chosen.  Hieftje  and  Bastiaans  (10)  have  shown  experimentally, 
through  measurement  of  the  rate  of  analyte  appearance  from  indiviu.il  solu- 
tion droplets  in  laminar  analytical  flames,  that  of  radius  of  a vaporizing 
solute  particle  follows  a relationship  of  form 

t * it  “ rd>  for  0 S t S ^ ' (4) 

k 

where  rQ  is  the  initial  dry  aerosol  particle  radius  (Jim) , the  radius  at 
time  t^,  and  1<  is  the  vaporization  rate  Omr/s).  Obviously,  the  mass  lost 
by  the  condensed-phase  particle  as  its  radius  decreases  must  be  gained  by 
the  vapor  phase.  Because  a certain  fraction,  of  that  mass  entering  the 
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vapor  phase  will  be  in  the  form  of  free  analyte  atoms,  the  relation  Co  • 
the  flux  of  atoms  (atoms/s)  into  the  vapor  phase,  >l'(t),  can  he  derived: 

4>(t)  = 2irp£k(r^  - kt)1/2  for  0 i t <.  r^/k  ('.) 

$(t)  = 0 for  t > r^/k  ((») 

o 

3 

where  ja  is  the  density  of  analyte  atoms  (atoms/V'm  ) in  the  condensed  phase. 
The  shape  of  this  vaporization  function  is  graphically  presented  as  curve 
A of  Figure  1. 

By  substituting  the  observed  analyte  vaporization  rate,  equation  5, 
into  the  solution  of  Fick's  Law,  (equation  3)  an  equation  can  be  derived 
which  describes  the  atom  concentration  as  a function  of  radial  distance 
in  the  vapor  cloud  and  of  time  from  the  initiation  of  vaporization: 


C(r,t) 


(6) 


Unfortunately,  equation  6 cannot  be  solved  in  closed  form  but  must  be 
evaluated  numerically. 

Because  of  the  potential  utility  of  an  analytical  solution  to  this 
convolution  integral  (equation  3),  it  is  worthwhile  to  npproxiatc  the 
vaporization  function  $(t)  by  simpler  but  similar  functions.  Two  such 
approximations  will  be  considered  here.  The  first  and  simplest  approxima- 
tion to  the  vaporization  function  is  one  representing  a constant  rate  of 
mass  loss;  i.e.,  it  assumes  that  analyte  atoms  are  liberated  at.  a constant 
rate  over  the  entire  vaporization  period:  < 

$(t)  = kv  for  0 s t s tv 


(7) 
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The  solution  of  equation  3 using  this  constant  vaporization  model  is  as 
follows : 


c(r't}  = EffiF  erfc  (sip) 


for  t £ t. 


and 


(8) 


c(r’t)  * 5S5F  erfc  (sm^) ' erf0  (a[b(t-w))>«) 

for  t > ty 

A closer  approximation  to  the  experimentally  observed  analyte  vapori- 
zation function  is  one  in  which  the  rate  of  vaporization  decreases  linearly 
with  time  during  the  entire  time  over  which  vaporization  takes  place.  The 
mathematical  statement  of  this  model  is  as  follows: 


#{t)  = qt  + Q 


for  0 s t s 


(9) 


where  £ and  £ are  the  slope  and  intercept  of  the  linear  function.  A graphical 

I 

comparison  of  the  true  vaporization  function,  equation  5,  and  these  approxi- 
mations, equations  7 and  9,  is  given  in  Figure  1.  The  solution  of  the  con- 
volution integral,  equation  3,  with  this  linear  approximation  to  the  vapori- 
zation rate,  equation  9,  is: 


C^t)  = Im?D 


for  0 s t s tj 


and 


(30) 


C(r't}  = & 


t + s + 

q 2D 


erfC(2lD?7V;?)-  6rfC  (2[D(t-tf)]^) 


TtJdT^ 


(t-tf)V2  e-r2Al>(t-tf)  _ tVSe-r2/4Dt 


for  t > tf 
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RESULTS  AND  DISCUSSION 
a/VWWWVV/WWWWWVV 

The  equations  developed  in  the  proceeding  section,  describing  the 
spatial  distribution  of  atoms  surrounding  a volatilizing  solute  particle, 
are  sufficiently  complicated  that  an  understanding  of  their  physical  rel- 
evance cannot  be  gleaned  by  simple  inspection.  Moreover,  a simple  graph  of 
radial  concentration  of  analyte  at  any  particular  time  is  difficult  to 
interpret.  Discerning  from  such  a plot  the  total  number  of  atoms  in  a 
particular  vapor  cloud  is  not  simple;  one  must  integrate  the  product  of  the 
graph's  value  times  a factor  representing  the  volume  of  a standard  (fixed 
thickness)  spherical  shell  at  each  radius.  Another  problem  with  this  type 
of  graph  is  that  the  model  employed  here  requires  that  the  concentration 
of  vapor  phase  analyte  be  infinite  at  the  center  of  the  vapor  cloud  during 
analyte  volatilization.  This  problem  derives  from  the  fact  that  the  vapor 

i 

I 

phase  atoms  are  being  injected  into  the  system  in  an  infinitely  small  volume. 

A better  appreciation  for  the  mathematical  model  can  be  gained  by 
plotting  the  number  of  atoms  contained  in  a thin  spherical  shell  of  uniform 
thickness  as  a function  of  the  shell's  radial  distance  from  the  origin. 

Let  the  thickness  of  the  shell  be  Ar,  and  let  the  shell  extend  in  radial 
distance  from  the  center  of  the  vapor  cloud  between  r - Ar/2  and  r + Ar/2. 

The  volume  of  this  spherical  shell  (AV)  can  be  expressed  as 

AV  = 4nr2Ar  + ir(Ar)3/3  (11) 

The  relationship  between  concentration  and  the  number  of  atoms  in  the  thin 
shell  now  becomes 


N(r,t)  = AV  C(r , t) 


(12) 


The  parameter  Ar  will  be  fixed  for  any  given  representation  and  must  be 


chosen  on  the  basis  of  the  resolution  desired.  A plot  of  N(r,t)  versus 

radial  distance  in  the  vapor  cloud  can  now  easily  be  UBed  to  determine  the 

number  pf  atoms  in  a sphere  of  a certain  sire,  or  Indeed  in  the  entire 

vapor  cloud,  simply  by  observing  the  area  under  the  curve  from  the  origin 

to  the  radius  desired.  A plot  of  equation  12  reveals  that  N(r,t)  approaches 

zero  as  the  radius  of  the  shell  becomes  small.  This  behavior  does  not  Imply 

that  concentrations  are  not  great  near  the  origin?  rather,  the  volume  of 

the  shells  near  the  origin  are  simply  so  small  that  vety  few  total  atoms 

reside  in  each,  even  at  high  atom  concentrations. 

To  evaluate  the  utility  of  the  models  described  in  the  previous  section, 

let  us  theoretically  examine  the  spatial  distribution  of  atoms  predicted 

from  s typical  aerosol  droplet  used  in  a practical  analytical  situation. 

Let  this  average  droplet  be  10  Vm  in  diameter  and  be  comprised  of  an  aqueous 

solution  containing  100  bg/ml  of  Ca  as  the  chloride  and  be  vaporised  in  on 

air/C^H^  flame.  When  desolvotion  of  this  droplet  is  complete  it  can  be 

assumed  that  the  resulting  aerosol  particle  will  he  a sphere  of  CnG^  of 

10  3 

radius  rQ  ■ 0.253  bm  and  density  p « 1.167  x 10  Ca  atoms/bm  (14)  and 
which  travels  essentially  at  the  flame's  velocity.  When  this  particle 
vaporizes,  it  is  expected  to  produce  a cloud  of  atoms  whose  spatial  charac- 
teristics change  with  time  in  accordance  with  equations  6 and  12.  However, 
two  parameters,  the  vaporisation  constant  and  the  diffusion  coefficient, 
need  to  be  defined  in  order  to  evaluate  these  equations.  Bastiaans  and 
Hieftje  (12)  have  shown  experimentally  that  Ca  vaporizes  in  an  air/^Hj 

flame  according  to  equation  .5  and  that  the  vaporization  constant  Jt  is 
2 -1 

1148  ym  s , Sue 11 man  (13)  has  measured  the  diffusion  coefficient  of  Ca 

a ? • 1 

in  an  air/CgH,,  flame  and  found  P - 3.0  x 10  urns  . Figure  2 allows  a 
three-dimensional  plot  of  the  predicted  atomic  population  per  alu*ll,  N ( r.O, 
versus  the  shell's  radial  distance  from  the  origin  of  the  vapor  cloud  am! 
vermin  time  for  this  typical  droplet  lu  an  nlr/OJl,,  flame. 
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In  calculating  Figure  2,  it  is  assumed  that  height  in  the  flame  is 
proportional  to  time  because  the  aerosol  particle  and  its  volatilization 
products  move  at  the  flame's  velocity.  At  a typical  rise  velocity  for 
an  air/C2H2  flame  of  10  m/s  (12,16,17),  one  millisecond  of  elapsed  time 
represents  an  upward  translation  of  one  centimeter  for  the  center  of  the 
vapor  cloud.  Consequently,  for  the  particle  under  consideration,  vapori- 
zation occurs  during  the  first  55.8  ys,  during  which  time  the  cloud  and 
particle  move  upward  0.558  mm.  Figure  2 shows  that  during  this  vapori- 
zation time,  the  number  of  liberated  analyte  atoms  increases  while  the 
diffusive  forces  spread  the  vapor  cloud  radially.  After  the  completion  of 
vaporization,  the  number  of  atoms  in  the  cloud  remains  constant  but  diffu- 
sion continues  to  spread  the  vapor  cloud. 

The  results  portrayed  in  Figure  2 are  not  unexpected.  However,  of 
greater  interest  is  what  happens  to  the  analyte  vapor  cloud  when  a matrix 
component  causes  a change  in  the  analyte  volatilization  rate.  In  the  mathe- 
matical models  proposed  here,  vaporization  is  an  important  factor  in  deter- 
mining the  spatial  distribution  of  analyte  atoms.  To  underscore  this  fact, 

Figure  3 depicts  the  differences  in  the  time-varying  atomic  spatial  distri- 
bution which  would  result  if  only  the  vaporization  rate  of  the  particles 

ft 

should  change.  Figures  3a,  b,  c,  and  d represent  radial  atomic  populations 
for  10  ys,  50  ys,  70  ys,  and  200  ys  after  the  initiation  of  analyte  volati- 

! 

lization  respectively.  In  each  instance  curve  A represents  the  number  of 
atoms  per  shell  (shell  width,  Ar,  is  100  ym  for  the  same  "typical"  CaC^ 

solution  droplet  used  to  calculate  Figure  2.  The  curves  labeled  £ and  £ j 

represent  similar  calculations  for  the  same  solution  droplet  but  with  a 
different  vaporization  rate  being  assumed.  Specifically,  curves  £ and  £ 
represent  a 20%  change  (slower  and  faster  respectively)  in  the  analyte  vapori- 

j 

t 


zation  rate. 
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In  Figure  3a,  one  sees  that  the  major  initial  effect  of  clung  lug  the 
vaporization  rate  constant  is  that  the  total  number  of  atoms  present  In  the 
vapor  phase  is  altered.  This  conclusion  derives  from  the  different  areas 
but  similar  shapes  of  the  three  curves  in  Figure  3a.  For  example,  comparing 
the  areas  of  curves  A and  B in  Figure  3a,  one  calculates  for  this  particular 
system  that  a particle  with  a 20%  slower  vaporization  rate  (curve  B)  yields 
about  20%  fewer  freed  analyte  atoms  after  10  ys.  In  comparison,  Figure  3b 
shows  similar  differences  in  the  number  of  atoms  for  the  three  vaporization 
rates  near  the  end  of  the  vaporization  process.  By  this  time.,  50  ys  after 
the  start  of  analyte  volatilization,  subtle  differences  in  the  shapes  of  the 
vapor  clouds  begin  to  appear.  The  apparent  trend  is  for  the  more  slowly 
vaporized  analyte  to  form  a more  compact  vapor  cloud.  This  trend  is  more 
dramatically  shown  in  Figure  3c,  where  analyte  volatilization  is  complete 
even  for  the  slowest  case  being  considered.  In  Figure  3c,  curves  A,  B,  and 
C have  equal  areas,  indicating  that  there  are  equal  numbers  of  atoms  in  the 
vapor  phase  after  volatilization  is  complete.  On  the  other  hand,  curve  B 
is  considerably  shifted  toward  the  origin,  demonstrating  that  slowed  vapori- 
zation of  analyte  causes  the  resulting  vapor  cloud  to  be  less  diffuse. 

Figure  3d  illustrates  that  at  a relatively  long  time  after  the  completion 
of  vaporization,  the  effects  caused  by  small  changes  (±  20%)  in  the  analyte 
vaporization  rate  are  still  present  although  of  lesser  magnitude  than  earlier 
in  the  cloud's  history. 

In  contrast  to  the  changes  wrought  by  an  altered  vaporization  rate 
(Figure  3),  Figure  4 reveals  the  effects  of  changing  analyte  diffusion 
coefficient  on  its  distribution  in  the  vapor  cloud.  The  curves  labeled  A 
in  Figure  4 are  again  the  predicted  time-varying  atomic  spatial  distri- 
butions for  the  same  "typical"  CaC^  solution  droplet  used  in  Figures  2 
and  3.  -In  all  cases,  curves  labeled  D and  1!  represent  a 20%  decrease  and 
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Increase,  respectively,  in  the  diffusion  coefficient  from  the  standard  C.iC I 
solution  droplet  plotted  in  curves  A.  Figures  4a,  b,  c,  and  d represent  the 
distribution  of  the  vapor-phase  Ca  after  10  ps,  50  ps,  /()  ps , and  ;’(H)  ps, 
respectively,  from  the  start  of  analyte  vaporizat  ion.  in  contrast  i <>  Figure 
3,  at  every  time  being  considered  here,  curves  A,  D,  and  K have  equal  area.  , 
indicating  that  the  vapor  clouds  being  represented  contain  an  equal  number 
of  atoms.  As  expected,  at  each  selected  time,  the  general  shape  of  curves 
A,  D,  and  E arc  similar.  Curve  D,  with  its  more  mobile  atoms,  consistently 
displays  a more  diffuse  vapor  cloud. 

Thus,  it  appears  that  both  vaporization  rate  and  diffusion  coefficient 

are  important  factors  in  determining  the  spatial  distribution  of  analyte 

atoms  in  a flame.  However,  the  effects  of  the  two  factors  are  distinguish- 
» 

able.  One  difference  in  the  effects  is  apparent  during  the  analyte  vapori- 
zation period,  i.e.,  before  the  time  when  vaporization  is  complete.  Figures 
3a  and  3b  show  that  when  vaporization  rates  differ,  the  total  number  of  atoms 
present  at  any  one  time  during  vaporization  will  also  differ.  On  the  other 
hand.  Figures  4a  and  4b  show  that  if  diffusion  coefficients  change,  only  the 
spread  of  atoms  and  not  their  number  will  change  at  a particular  time  during 
volatilization. 

Figures  3d  and  4d  show  that  at  long  times  after  analyte  volatilization 
is  complete,  another  difference  exists  between  the  effects  of  changing 
vaporization  rates  and  changing  diffusion  coefficients.  The  variation  in 
spatial  spread  of  atoms  caused  by  alterations  in  analyte  volatilization  rate 
become  relatively  small,  as  seen  in  Figure  3d.  In  contrast,  the  dependence 
of  atomic  distribution  on  vapor-phase  analyte  mobilities  can  be  seen  in 
Figure  4d  to  be  still  relatively  large.  At  intermediate  times  the  effects 
of  the  two  changes  are  very  similar,  as  can  be  seen  by  comparison  of 
Figures  3c  and  4c. 


Because  the  spatial  distribution  of  atoms  is  controlled  by  both  tin- 
diffusion  coefficient  and  the  vaporization  rate,  the  mechanism  for  any 
"lateral  diffusion  interference"  might  result  from  changes  in  either  or 
both.  However,  in  conventional  flame  systems,  it  will  be  difficult  to 
ascribe  the  observed  interference  to  either-  source.  While  differences  in 
these  effects  seem  dramatic  in  Figures  3 and  4,  the  unpredictability  and 
lack  of  control  of  typical  nebulizers  and  burners  makes  it  unlikely  that 
the  effects  of  diffusion  and  vaporization  will  be  distinguishable  from  each 
other.  The  differences  noted  during  vaporization  would  require  optical 
resolution  on  the  order  of  hundreths  of  millimeters  in  a conventional 
flame  system.  Other  differences,  noted  in  Figures  3 and  4 long  after 
particle  volatilization  is  complete,  would  also  be  difficult  to  observe 
in  typical  analytical  flames.  Such  observations  would  have  to  be  made 
high  in  the  flame,  where  the  flame  is  unstable  and  where  atmospheric  entrain 
ment  becomes  significant. 

Convenient  incorporation  of  the  foregoing  theory  into  models  for  atom 
formation  and  into  studies  on  interference  effects  requires  a closed-fonn 
solution  to  the  integral  equations;  in  turn,  such  closed-form  solutions 
require  approximations  to  the  vaporization  function  $(t)  to  be  employed. 
Because  equations  8 and  10  represent  such  approximate  solutions,  let  us 
evaluate  how  well  they  represent  the  true  (equation  6,  numerically  inte- 
grated) behavior  under  a variety  of  conditions  and  which  would  be  best  to 
employ  in  further  treatments.  In  the  constant  vaporization  (CV)  model, 
the  vaporization  rate,  kv>  and  vaporization  period,  ty,  can  be  varied  to 
produce  the  best  fit  to  the  true  vaporization  characteristics  of  the  analyte 
represented  by  equation  3 and  portrayed  in  curve  A of  Figure  1.  The  best 

fit  is  found  by  keeping  t equal  to  the  true  function's  vnporiz.it  Ion  period, 
2 

r D/k.  Requiring  that  both  vaporization  functions  produce  the  same  number 
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of  atoms  sets  the  value  of  the  vaporization  rate,  ky.  In  the  case  of  the 
10  Jim  droplet  of  CaC^  used  in  Figures  2 through  4,  a best  fit  of  the  CV 
model  to  the  true,  experimentally  confirmed  vaporization  function  is: 


(13) 

In  the  linear  vaporization  model  (LV) , three  parameters  are  available 
to  adjust  the  approximate  vaporization  function:  q,  Q,  and  t^.  The  values 
of  these  parameters  are  determined  by  minimizing  the  integral  of  the  square 
of  the  difference  between  the  true  and  approximate  functions.  This  model 
best  fits  the  experimentally  determined  vaporization  of  the  10  pm  droplet 
of  CaCl^  when: 

(14) 

i 

I 

Figure  1 graphically  compares  the  true,  experimentally  confirmed  vaporization 
function  with  these  two  approximation::,  represented  by  equations  13  and  14. 

Using  the  above  approximations  to  the  vaporization  function.  Figure  5 
compares  the  constant  vaporization  (CV)  and  the  linear  vaporization  (LV) 
models  to  the  behavior  predicted  by  the  more  rigorous  theory  of  equation  6. 
Figure  5a  portrays  the  atomic  spatial  distributions  predicted  for  10  ys  after 
the  initiation  of  analyte  volatilization.  Both  LV  and  CF  approximations  pre- 
dict the  correct  shape  for  the  atomic  vapor  cloud  at  this  very  early  time, 
even  though  the  CV  model  fails  to  predict  the  correct  number  of  atoms  in  the 
vapor  cloud.  Figure  5b  shows  that  near  the  end  of  vaporization  the  constant 
vaporization  model  more  closely  predicts  the  correct  number  of  atoms  in  the 
vapor  cloud  but  the  atoms  are  decidedly  closer  to  the  origin  than  the  true 
function  indicates;  at  this  same  time,  the  linear  vaporization  function  fits 
the  true  behavior  very  well  indeed.  Figure  5c.  shown  that  shortly  after  the 
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completion  of  vaporization  the  linear  vaporization  inode]  Is  indistinguishable 
from  the  numerical  integration  of  equation  6 while  the  constant  vaporization 
model  predicts  a considerably  smaller  vapor  cloud.  A relatively  long  time 
after  the  completion  of  vaporization,  portrayed  in  Figure  5b,  both  the 
constant  and  linear  vaporization  models  provide  a reasonable  approximation 
to  the  true  predicted  behavior.  Thus  it  appears  that  the  linear  vaporization 
model,  represented  by  equation  10,  can  be  used  to  reasonably  approximate 
atomic  distribution  in  e large  number  of  instances  while  the  constant  vapori- 
zation model  only  approximates  the  spatial  atomic  distributions  at  long  times 
after  vaporization  is  complete. 

CONCLUSIONS 

'WVWWVWV 

• The  spatial  distribution  of  atoms  created  by  the  vaporization  of  a single 
aerosol  particle  in  a laminar  flame  can  be  mathematically  modeled.  Equation 
3 is  the  general  mathematical  representation  of  this  model.  A rigorous  solu- 
tion of  the  model,  using  an  experimentally  determined  vaporization  rate, 
requires  numerical  integration  of  equation  6.  The  results  of  this  integra- 
tion have  been  used  to  prove  that  the  rate  of  volatilization  does  indeed 
play  a significant  role  in  determining  the  spatial  distribution  of  atoms  in 
an  individual  vapor  cloud. 

The  theory  of  spatial  distribution  of  analyte  atoms  vaporized  from 
individual  aerosol  droplets,  as  presented  here,  distinguishes  between  the 
effects  of  analyte  mobility  and  vaporization.  For  this  reason,  these  theories 
should  be  useful  in  determining  the  mechanism  of  the  "lateral  diffusion" 
interferences  in  practical  flame  atomic  absorption  spectrometry. 

Coupled  with  a prediction  of  the  distribution  of  aerosol  droplets 
entering  a practical  flame,  this  theory  should  aid  in  predicting  the  spatial 
distribution  of  atoms  in  even  complex  flame  systems  (18) . Such  accurate 
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descriptions  of  analyte  distributions  will  be  essential  in  the  quest  to 
achieve  absolute  flame  atomic  absorption  measurements. 

In  some  applications  of  this  theory,  an  equation  which  does  not  require 
numerical  integration  may  be  desirable.  Such  analytical  solutions  are  pro- 
vided by  equations  8 and  10  through  the  use  of  approximations  to  the  vapori- 
zation rate.  While  these  solutions  are  not  as  rigorous  as  equation  6,  the 
error  they  cause  under  various  conditions  has  been  evaluated.  The  linear 
vaporization  model,  equation  10,  can  be  used  to  reasonably  approximate  atomic 
distributions  both  during  and  after  analyte  vaporization.  While  the  constant 
vaporization  model,  equation  8,  never  predicts  spatial  distributions  as 
accurately  as  the  linear  model,  the  errors  become  small  at  long  times  after 
vaporization.  For  some  applications,  the  relative  mathematical  simplicity 
of  equation  8 may  be  more  important  than  the  errors  that  its  use  will  cause. 

Of  course,  the  calculation  of  spatial  atomic  distributions  for  any 
particular  system  requires  a knowledge  of  the  analyte  diffusion  coefficient 
and  vaporization  rate.  Bastiaans  and  Hieftje  (12)  and  Clampitt  and  Hieftje 
(19)  have  experimentally  studied  rates  of  analyte  vaporization  in  flames. 

The  diffusion  rates  of  atoms  in  flames  have  been  measured  by  several 
investigators  (15,20). 

The  theory,  as  presented  here,  is  limited  in  that  it  does  not  account 
for  perturbations  in  the  atomic  vapor  cloud  shape  due  to  the  overlap  of 
vapor  clouds  from  neighboring  aerosol  particles.  This  theory  also  does  not 
account  for  the  effects  of  chemical  reactions  which  may  change  the  free 
atomic  populations  in  different  portions  of  the  vapor  cloud. 

In  a subsequent  publication  (21),  this  theory  will  be  tested  with 
experimentally  measured  spatial  distributions  of  analyte  vaporized  from 
individual  solution  droplets  injected  Into  an  analytical  flame.  These  mea- 
surements will  be  used  to  deduce  effective  analyte  diffusion  coefficients 
and  investigate  the  lateral  diffusion  interference. 
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GLOSSARY 
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C(r,t)  concentration  of  analyte  atoms  in  the  vapor  phase 
CV  constant  vaporization  model 

k vaporization  rate  constant  for  analyte  atoms 

kv  vaporization  rate  for  the  constant  vaporization  model 

LV  linear  vaporization  model 

N number  of  vapor  phase  analyte  atoms  in  a spherical  shell  of 

thickness  Ar 

Nt  total  number  of  vapor  phase  analyte  atoms  in  a single  vapor  cloud, 

q slope  of  linear  vaporization  model 

Q intercept  of  linear  vaporization  model 

r distance,  from  the  center  of  the  atomic  vapor  cloud 

r^  time  dependent  radius  of  the  vaporizing  analyte  particle 

rQ  initial  radius  of  the  completely  desolvated  aerosol  particle  I 

Ar  thickness  of  a spherical  shell  within  an  atomic  vapor  cloud 

t time,  measured  from  the  onset  of  analyte  vaporization 

tj  integration  limit  for  convolution 

t time  for  completion  of  particle  vaporization 

t’  delay  parameter  used  in  the  convolution  integral 

Av  volume  of  a spherical  shell  of  atomic  vapor 

$(t)  vaporization  function  (atoms/s)  for  analyte  atoms 

3 

p density  of  analyte  atoms  (atoms/pm  ) in  the  condensed  phase 


fraction  of  the  vapor  phase  analyte  atoms  which  are  neither  bonded 
to  other  atoms  or  ionized. 
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FICURi;  CAPTIONS 

WjWi/WI/WWVI 


Figure  1.  Experimentally  verified  and  approximate  models  for  volatilization 
of  CaCl2  In  an  air-acctylcne  flame.  Curve  A allows  the  rate  of 

volatilization  of  Ca  atoms  from  a sphere  of  CaCl„  density 

10  3 ^ 

p = 1.16  x 10  Ca  atoms/pm  , liaving  an  initial  dry  radius 

r = 0.253  pm  and  vaporizing  in  an  air-acetylene  flame  accord- 

® 2—1 
ing  to  equation  5.  Vaporization  constant  is  k = 1148  s 

Reference  point  of  the  time  axis,  t = 0,  is  the  initiation  of 
Ca  volatilization.  Curve  CV  is  the  flux  of  Ca  atoms  into  the 
vapor  phase,  predicted  for  the  same  physical  system  by  the  con- 
stant vaporization  model,  equation  13.  Curve  LV  is  the  approxi- 
mation to  the  experimentally  confirmed  volatilization  of  Ca 
atoms,  curve  A,  by  the  linear  vaporization  model,  equation  14. 

Figure  2.  Radial  and  time  dependent  vapor-phase  Ca  atom  populations.  The 
population  of  vapor-phase  Ca  atoms  in  spherical  shells  of  uni- 
form thickness,  Ar  = 100  pm  is  plotted  with  respect  to  the  dis- 
tance from  the  center  of  the  vapor  cloud  to  the  spherical  shell 
and  with  respect  to  time.  The  time  reference  point,  t = o,  is 
set  at  the  onset  of  Ca  atom  vaporization.  The  hypothetical 
aerosol  droplet  which  produced  this  distribution  of  vapor-phase 
Ca  atoms  was  a 10  pm  diameter  droplet  containing  100  pg/ml  Ca 

as  CaC^.  After  desolvation,  vaporization  of  Ca  was  calculated 

2 -1 

according  to  equation  4 with  k = 1148  pm  s and  is  shown 

graphically  in  Figure  1.  Ca  vaporization  is  complete  after 

first  55.8  ps  plotted.  The  synergic  effects  of  vaporization 

2 

rate  and  Fick's  Law  diffusion,  D = 3.0  cm  /s,  in  creating  the 
spatial  distribution  of  atoms  shown  above  is  given  in  mathe- 
matical form  in  equation  6. 


Figure  3.  Changes  in  analyte  spatial  distribution  caused  by  changes  in 
particle  volatilization  rate.  Curves  labeled  A^  represent  the 
Ca  atoms  in  a 100  pm  thickness  spherical  shells  from  the  same 
aerosol  droplet  assumed  for  Figures  1 and  2.  Radial  populations 
are  predicted  by  the  combination  of  equations  6 and  12.  Curves 
labeled  _B  and  are  Ca  atom  populations  produced  from  identical 
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Figure  4. 


Figure  5. 


droplets,  but  the  droplet  in  H vaporizes  202  more  slowly,  and 
£ vaporizes  202  more  rapidly  than  the  droplet  in  A.  The  effects 
of  vaporization  rate  changes  vary  with  time.  Pl.cs  a,  b,  c, 
and  d of  this  figure  represent,  respectively,  10  ps,  50  ps, 

L» 

70  ps,  and  200  ps  after  the  initiation  of  Ca  volatilization. 

Note  differences  in  both  horizontal  and  vertical  scale  in 
parts  a,  b,  c,  and  d. 

Changes  in  analyte  spatial  distribution  caused  by  changes  in 
its  diffusion  coefficient.  The  curves  labeled  A represent  the 
Ca  atoms  in  100  pm  thick  spherical  shells  from  the  same 
aerosol  droplet  used  in  Figures  1 and  2 and  the  curves  labeled 
A in  Figure  3.  The  curves  labeled  D and  E represent  system 
identical  to  that  used  for  curve  A with  only  the  increase  and 
decrease,  respectively,  in  the  diffusion  coefficient.  The 
relative  magnitude  of  the  change  in  the  diffusion  coefficient 
is  identical  to  the  relative  magnitude  of  the  change  in  vapori- 
zation rate  used  in  Figure  3 (i.e.,  20%).  Parts  a,  b,  c,  and 
d represent  the  radial  analyte  population  at  10  ps,  50  ps,  70 
ps,  and  200  ps  after  the  start  of  Ca  vaporization.  Note  scale 
changes  in  different  figures. 

Comparison  of  the  predicted  radial  analyte  population  produced 
by  the  use  of  approximate  vaporization  functions.  The  curves 
labeled  A represent  calculations  identical  to  those  of  Figure  2 
and  curves  A of  Figures  3 and  4.  The  curves  labeled  CV  repre- 
sent the  constant  vaporization  approximation  to  this  model 
(cf.  equations  8,  12,  and  13).  The  curves  labeled  LV  repre- 
sent the  linear  vaporization  approximation  to  this  model  used 
for  the  curves  A (cf.  equations  10,  12,  and  14).  The  differ- 
ences in  the  vaporization  rates  of  these  three  models  are  dis- 
played graphically  in  Figure  1. 
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